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Abstract 


Recent  experiments  have  disclosed  electrical  effects  well  ahead  of  an 
advancing  shock  front.  Certain  of  these  effects  have  been  attributed  to  the 
diffusion  of  electrons  through  the  shock  front  from  the  ionized  region  behind 
it.  In  this  report,  the  diffusion  hypothesis  is  examined  in  terms  of  a  simple 
heuristic  model,  in  which  the  diffusion  is  assumed  to  take  place  from  a  plane 
electron  source  moving  with  the  shock  velocity.  The  predictions  of  this  model 
are  shown  to  be  generally  consistent  with  the  available  experimental  evidence. 
It  is  further  shown  that  the  transient  behavior  of  the  precursor  electron  dis¬ 
tribution  is  capable  of  increasing  the  Doppler  shift  of  an  electromagnetic 
wave  reflected  from  the  shock.  Although  the  transient  decays  rapidly  after 
the  creation  of  the  shock,  there  are  conditions  under  which  it  should  last 
long  enough  to  allow  an  experimental  check  of  this  effect. 
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I.  Introduction 


Recent  experiments  in  shock  tubes  have  disclosed  the  presence  of  a 
variety  of  electrical  effects  in  the  region  well  ahead  of  the  advancing 
shock  front. ^  ^  Several  hypotheses  have  been  proposed  to  explain  the 
existence  of  these  effects,  but  no  one  of  them  alone  can  account  for  all 
of  the  observations.  The  observations  themselves,  which  consist  of  sig¬ 
nals  recorded  from  probes  placed  along  the  shock  tube,  may  be  classified 
according  to  the  apparent  velocity  of  the  precursor  signal.  Soane  of  these 
signals  appear  to  propagate  with  effectively  infinite  velocities,  and  are 
presumably  caused  by  photo-effects  from  the  delayed  luminosity.  ^  Others 
propagate  with  velocities  of  the  order  of  the  shock  speed,  and  have  been 
attributed  to  the  diffusion  of  electrons  upstream  through  the  shock  front. ^ 

Signals  in  the  latter  class  have  been  observed  even  when  the  delayed  lumi- 

(2) 

nosit y  did  not  appear''  ,  and  seem  to  have  been  detected  only  when  the  shock 

was  formed  in  one  of  the  noble  gases.  It  is  significant  that  the  noble 

gases  exhibit  a  pronounced  Ramsauer  effect  for  electron  energies  corresponding 

to  the  temperatures  generally  existing  in  the  equilibrium  region  behind  a 

hypersonic  shock  front.  The  recent  discovery  of  transient  "double  fronts" 

fin 

in  small  bore  electric  shock  tubes'^'  might  be  interpreted  as  further  evi¬ 
dence  of  precursor  electron  activity,  although  this  behavior  has  been 
ascribed  to  the  high-energy  electrons  from  the  energizing  electric  arc^^, 
rather  than  to  the  ionization  created  by  the  shock. 

At  present  the  diffusion  hypothesis  rests  almost  entirely  on  somewhat 
meager  experimental  evidence,  the  complex  character  of  the  shock  front  region 
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making  an  adequate  theoretical  investigation  extremely  difficult.  We  will 
not  attempt  such  an  investigation  here.  Instead,  we  will  review  the  current 
picture  of  how  ionization  takes  place,  then  speculate  briefly  on  some  aspects 
of  the  diffusion  problem.  In  so  doing,  it  is  hoped  that  the  simple  heuristic 
model  on  which  the  conclusions  of  this  report  are  based  will  be  made  plau¬ 
sible  . 

It  is  reasonable  to  ask  first  about  the  electron  and  ion  densities  which 
would  appear  if  diffusion  were  ignored.  According  to  Petschek  and  Byron v ' ' 
(among  others),  the  ionization  process  can  be  roughly  divided  into  two  stages. 
The  shock-heated  gas  reaches  thermal  equilibrium  directly  behind  the  shock 
front,  and  as  one  moves  downstream,  a  certain  amount  of  initial  ionization 
takes  place  *8  a  result  of  atom-atom  and  atom- impurity  collisions.  As  the 
electron  density  grows,  electron-atom  collisions  take  over  as  the  primary 
ionizing  mechanism,  and  the  gas  is  fairly  quickly  ionized  to  the  equilibrium 
level.  Thus  we  have  the  picture  of  the  electron  (and  ion)  concentration 
starting  from  zero  at  the  shock  front,  rising  slowly  behind  the  shock  front 
for  a  while,  then  increasing  rapidly  to  the  equilibrium  level.  Although  the 
details  will  depend  on  the  type  of  gas,  the  shock  speed,  the  impurity  con¬ 
centration,  etc.,  the  general  behavior  of  the  plasma  concentration  may  be 
represented  by  the  following  picture: 
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This  does  not,  of  course,  represent  the  steady-state  electron  concentration, 
since  diffusion  has  been  ignored.  The  electrons  have  much  larger  diffusion 
coefficients  than  the  ions,  so  they  will  tend  to  diffuse  out  of  the  plasma, 
leaving  the  ions  behind.  The  electrostatic  space-charge  fields  arising  from 
this  charge  separation  will  retard  the  flow  of  electrons,  but  we  cannot 
predict  to  what  extent  and  in  what  manner  without  first  producing  a  convincing 
solution  for  the  diffusion  problem  associated  with  the  ioniaation  variations 
on  which  Fig.  1  is  based.  It  is  just  here  that  the  theory  is  stalled.  Electron 
flow  in  the  presence  of  space-charge  effects  has  been  considered  in  other  con¬ 
texts,  however,  and  it  will  be  helpful  to  speculate  from  the  point  of  view  of 
a  few  of  these  results. 

We  observe,  first  of  all,  that  the  zeroth-order  concentrations  change 
from  very  small  values  directly  behind  the  shock  front  to  very  large  values 
in  the  region  of  ionization  equilibrium  behind  the  "onset  point".  Now  the 
two  extremes  of  electron  diffusion  out  of  a  plasma  concentration  may  be 
characterized  by  the  ambipolar  and  the  free  diffusion  coefficients,  correspond¬ 
ing  respectively  to  high  and  low  plasma  densities.  Where  there  is  a  large 
variation  of  plasma  density,  as  in  the  case  being  considered  here,  one  may 
expect  a  transition  between  the  ambipolar  and  free  diffusion  regimes,  and 
this  transition  may  affect  the  character  of  the  entire  diffusion  process  in 
a  complicated  way.  For  example,  Allis  and  Rose^  have  considered  transition 
diffusion  in  the  case  of  a  microwave  discharge  in  a  cavity,  and  have  obtained 
good  agreement  between  theory  and  experiment.  Although  diffusion  into  an 
open  region,  as  through  a  shock  front,  introduces  additional  complications, 
the  above  suggests  that  the  diffusion  of  electrons  from  the  low-density  tail 
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of  the  distribution  shown  in  Fig„  1  might  be  governed  by  a  coefficient 
approaching  that  for  free  diffusion.  It  is  important  that  this  possibility 
should  exist,  because  if  electrons  are  to  succeed  in  outrunning  the  shock 
to  the  extent  suggested  by  the  experimental  data,  there  must  be  a  reservoir 
of  electrons  with  large  diffusion  coefficients  available  behind  the  shock 
front.  An  objection  can  be  raised  here  that  the  electrons  created  by  the 
ionization  processes  in  the  onset  region  are  left  with  relatively  low  energies, 
and  they  lack  sufficient  time  to  come  into  equilibrium  with  the  shock-heated 
gas.  Therefore  even  the  free  diffusion  coefficient  will  be  too  small  to 
account  for  the  observed  precursor  distributions. 

We  may  refer  here  to  an  analogous  situation  encountered  in  the  low-volt¬ 
age  arc.  In  this  device  it  is  found  that  a  stable  arc  in  a  noble  gas  can  be 
maintained  at  voltages  which  are  less  than  half  the  excitation  potential  for 
the  gas  involved.  Probe  measurements  disclose  the  existence  of  a  positive 
ion  concentration  in  the  center  of  the  arc  which  results  in  an  accelerating 
field  from  cathode  to  space-charge,  and  a  strong  decelerating  field  from 
space-charge  to  anode.  Electrons  from  the  cathode  are  therefore  accelerated 
to  the  excitation  potential  at  the  center  of  the  arc,  but  having  excited  the 
gas  atoms,  are  left  with  virtually  zero  energy.  The  magnitude  of  the  observed 
arc  current,  however,  requires  that  a  large  fraction  of  these  electrons  must 
overcome  a  potential  barrier  of  about  7v  to  reach  the  anode.  The  question: 

where  do  they  get  the  necessary  energy?  This  question  was  successfully  ans- 

(g) 

wered  by  Druyvesteyn,  '  who  showed  that  the  exhausted  electrons  gain  energy 
from  those  electrons  which  arrive  from  the  cathode,  but  .have  failed  to  excite 
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any  of  the  gas  atoms.  The  relaxation  time  for  this  equilibration  is  only 
about  one  mean  free  time  between  electron-electron  collisions,  which  is 
very  small  compared  with  the  large  times  required  for  a  low  energy  electron 
to  reach  equilibrium  with  the  hot  neutral  gas.  If  we  imagine  some  such 
process  to  be  operative  in  heating  the  cool  electrons  in  the  onset  region 
directly  behind  the  shock  front,  we  are  left  with  the  problem  of  finding  a 
source  of  hot  electrons  to  do  the  job.  This  is  not  hard,  because  the  only 
source  of  hot  electrons  lies  behind  the  onset  point,  in  the  region  of 
ionization  equilibrium.  The  supply  of  electrons  here  is  copious  and 
energetic,  so  there  now  exists  the  possibility  of  "thermionic  emission" 
from  this  region.  What  is  meant  by  this  is  that  electrons  from  the  high- 
energy  tail  of  the  Maxwell  distribution  move  across  the  potential  barrier, 
which  is  produced  by  charge  separation,  and  mingle  with  the  low  energy 
electrons  in  the  onset  region.  In  this  way  the  electrons  close  to  the 
shock  front  could  be  given  sufficient  energy  to  diffuse  upstream. 

The  considerations  of  the  last  two  paragraphs  suggest  a  simple  model 
from  which  certain  implications  of  electron  diffusion  may  be  found.  If  a 
"reservoir"  of  diffusible  electrons  is  to  exist  behind  the  shock  front,  then 
we  must  be  able  to  define  a  steady-state  density,  n^,  for  the  electrons  con¬ 
tained  in  it.  (The  reader  is  cautioned  not  to  take  the  term  "reservoir"  too 
literally.)  Since  the  probability  is  low,  on  a  random-walk  basis,  that  an 
electron  far  behind  the  front  will  contribute  to  the  precursor  distribution, 
we  may  assume  that  this  reservoir  is  located  in  a  narrow  region  of  width  d, 
directly  behind  the  shock  front.  In  order  to  maintain  the  steady-state 
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density  in  this  region,  the  shock  must  generate  new  electrons  at  the  rate 
n  vs/d  per  unit  volume  to  replace  those  that  are  left  behind  as  the  shock 
moves  along.  If  we  consider  a  one-dimensional  shock,  then  we  approximate 
this  picture  by  collapsing  the  precursor  reservoir  into  a  plane  electron 
source  o*f  strength  (n  v  / d)  x  (d)  =  n  v  ,  located  at  the  shock  front  itself. 

p  s  p  s 

We  will  denote  the  temperature  of  this  source  by  ,  In  this  way,  then,  we 

have  replaced  the  real,  but  intractable,  situation  by  a  fictitious  model 

consisting  of  a  plane  source  moving  with  the  shock  velocity  v  and  characterized 

s 

by  two  parameters,  n^  and  ,  into  which  are  lumped  the  unknowns  of  the  problem. 
If  accurate  measurements  could  be  made  of  the  precursor  distribution,  these 
parameters  could  be  determined  empirically. 

In  the  next  section  the  precursor  distribution  predicted  by  this  model 
will  be  derived.  Following  this,  it  will  be  shown  that  if  a  reasonable 
interpretation  of  the  available  measurements  is  allowed,  the  time-dependence 
of  the  predicted  precursor  behavior  is  consistent  with  the  experimental  results. 
Finally,  it  is  suggested  that  under  suitable  conditions,  an  electromagnetic  wave 
incident  upon  a  shock  front  might  experience  a  transient  enhancement  of  its 
Doppler  shift. 


II*  The  Precursor  Distribution 

Before  examining  the  transient  behavior  of  the  precursor  electron 
distribution,  let  us  obtain  the  steady-state  distribution  by  means  of  an 
elementary  argument.  Consider  a  shock-based  reference  system  in  which  the 
shock  front  is  fixed  and  the  driven  gas  streams  through  it  with  velocity 
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v  ,  as  in  Fig.  2  belows 
s 


Fig.  2 


Under  steady-state  conditions  in  this  system,  the  diffusing  electrons  must 
arrange  themselves  into  a  distribution  of  such  a  form  that  the  drift  velocity 
due  to  the  concentration  gradient  is  equal  and  opposite  to  the  velocity  of 
the  streaming  gas.  This  is  expressed  by  the  equation 


k~Tp  _l_  _  g- 

<rr\i/c  n  dS,  5 


(1) 


where  £  is  the  distance  upstream  measured  from  the  shock  front,  T  is  the 

P 

electron  temperature,  m  is  the  mass  of  the  electron  and  vq  is  the  collision 
frequency  for  electrons  of  temperature  in  the  neutral  gas.  (The  electro¬ 
static  self -field  of  the  distribution  has  been  ignored  in  view  of  the  rel¬ 
atively  low  electron  densities  expected.)  Since  it  is  being  assumed  that  these 


electrons  emerge  from  a  region  maintained  at  a  steady-state  density  n^  ,  Eq.  1 
may  be  solved  as 


=  ktpe 


(2) 


0  o 

where  D  =  v  /3v  (v  =  3kT  /m)  is  the  coefficient  for  free  diffusion  of 

O  C  O  p 

the  electrons  in  the  neutral  gas.  Weymann  claims  that  this  £ -dependence 
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can  be  verified  if  is  taken  as  the  equilibrium  temperature  behind  the 

shock  front. ^  if,  in  addition,  an  absolute  measurement  of  electron  density 

can  be  made  for  some  value  of  £  ,  then  n  can  be  determined.  This  is  a  model 

P 

parameter,  however,  and  has  no  particular  physical  significance. 

In  order  to  obtain  the  transient  behavior,  the  plane-source  model  de¬ 
scribed  at  the  end  of  the  last  section  will  be  used.  We  consider  a  one¬ 
dimensional  shock  which  is  created  at  time  t  =  0  at  the  origin  of  coordinates, 
x  ■  0.  Thus  one  might  have  in  mind  a  shock  tube  in  which  a  diaphragm  at  x  »  0 
is  burst  at  t  «  0.  We  will  be  interested  in  laboratory  measurements  of  pre¬ 
cursor  effects,  hence  in  the  electron  density  existing  at  a  fixed  point  x  at 
a  time  t  before  the  shock  front  has  arrived.  Figure  3  represents  a  "snap¬ 
shot"  of  the  advancing  shock  front,  showing  its  position  relative  to  the  point 
of  observation  at  the  time,  t,  of  observation. 


Fig.  3 

Observe,  however,  that  the  precursor  electrons  have  been  diffusing  out  of 
the  shock  front  over  the  entire  time  interval  from  0  to  t,  thus  the  density 
at  (x,t)  represents  an  accumulation  of  electrons  created  at  earlier  times. 
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This  accumulation  process  may  be  described  analytically  as  follows.  If  we 
let  x*  be  the  position  of  the  shock  front  at  a  time  t'<t,  (see  Fig.  3), 
then  in  the  time  interval  dt !  the  source  at  the  front  produces  n^vgdt’ 
potential  precursor  electrons.  Some  of  these  diffuse  toward  x,  and  their 
contribution  to  the  precursor  density  at  x  at  the  later  time  t  can  be 

written  - 

ftc-x') 


cWx,t;  X/i 0 


(3) 


where  D  is 
mulation  at 
from  0  to  t, 


the  diffusion  coefficient  defined  earlier.  The  total  accu- 

(x,t)  is  found  by  putting  x1  =  v  t1  and  integrating  Eq.  3 

s 


ft  Cx-Wtf 

^pvr?  (  *T>(*0  At' 

\l4-r-I)  \  . 


(U) 


The  integral  can  be  evaluated  by  converting  it  into  the  form  of  a  convolution 
of  two  functions.  In  terms  of  the  variable  t  *  t  -  t*  ,  Eq.  1*  becomes 


HCX/t) 


yyvTs 

MttV 


At 

'Tr*’ 


(5) 


Multiplying  out  the  exponent. 


-SI 

Cy-vrsfe-r)j 

^  Dr 


(x  -  'Xjt) 
4T>r 


and  using  the  notation  (x-vt)=£,we  may  write  {$)  as 

s 
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r\C<t)  - 


~*T>r 

i—  e41>  kr 
FF  S 


The  integral  in  (6)  is  now  in  the  form  of  a  convolution. 


4*%  =  ^  4 (?)<^(±-r)br  j 


which  may  be  expressed  as 


W  c*«. 


where  ^  and  ^  are  the  respective  Laplace  Transforms  of  f  and  g , 
and  ^  *  denotes  the  inverse  transform.  The  transforms  of  the  two  functions 
in  the  integrand  of  (6)  may  be  found  in  almost  any  table  of  Laplace  Trans¬ 


forms  to  be 


>  _ 


3- 


r  *1 T  -]  -i 

4  =  t[e<>  J  -  , 


,1  2 

where  q  -  (p/D)  2  and  a  =  v  /liD.  Thus,  if  we  let  I  be  the  integral  in 


Eq.  6,  we  may  write 


x  -  W 


.^Cp-sO 
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Performing  the  inversion,  one  obtains. 


which,  when  substituted  into  Eq.  6,  yields  n(x,t)  in  the  form: 


The  erfc  function  decays  rapidly  with  increasing  argument,  so  the  first  erfc 
in  the  brackets  of  (12)  is  always  much  larger  than  the  second  (except  in  a 
small  neighborhood  of  t  =  0.)  For  most  of  the  interesting  values  of  vg  and 
D,  the  first  term,  even  with  the  exponential  factor,  will  dominate  the  expres¬ 
sion,  so  we  will  approximate  n(x,t)  by 

«|nPe.  ”  (~j===r)  .  (13) 

It  is  of  interest  to  verify  that  this  expression  will  lead  to  the 
steady-state  distribution  obtained  earlier.  As  t  00  ,  the  argument  of 
the  erfc  goes  to  -  «  for  any  £  ,  and  the  erfc  approaches  2.  The  variable 
X  may  be  viewed  as  the  result  of  a  transformation  from  the  laboratory  system 
to  a  system  moving  with  the  shock,  so  Eq.  13  reduces  to  Eq.  2  in  the  limit  as 
t  becomes  infinite.  For  finite  values  of  t,  however,  we  see  that  the ’steady- 
state  distribution  is  "modulated"  by  the  time-dependent  erfc  function.  The 
form  of  n  given  by  Eq.  13  is  particularly  convenient  for  describing  the 
effect  of  this  modulating  factor,  since  it  can  be  considered  to  be  a 
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t-parametric  description  of  the  precursor  distribution  in  a  system  of  ref¬ 
erence  moving  with  the  shock. 

We  can  get  a  general  idea  of  what  to  expect  by  examining  the  extremes 

of  the  dependence  of  the  erfc  on  its  argument.  For  negative  values  of 

y<-l,  erfc  (y)  -V  2  very  rapidly,  while  for  positive  values  of  y>+l, 

_  2  _ 

erfc  (y)  ->  e~^  /y  y/rT  .  This  means  that  the  precursor  electron  density 

given  by  Eq.  13  decays  in  a  different  way,  according  as  £  >t  t,  That  is 
J  s 

iLd 

-  J>  ?  _ 

r\(it)  ->  r\P  a  /  i  <  vr*t  -  y  (Hi) 

(4 1  uy£ 

^  e  ^  /  i  >  vr5i  +  ' .  (15) 

The  transition  between  these  two  regimes  of  decay  may  be  taken  as  the  point 

for  which  £  »  vgt,  therefore  the  transition  point  is  moving  away  from  the 

shock  front  with  the  shock  velocity  v  .  In  the  laboratory  system  we  have 

s 

x  =  2v  t  ,  so  the  transition  point  moves  with  twice  the  shock  velocity, 
s 

One  is  tempted  to  regard  the  region  over  which  this  transition  takes  place 
as  a  kind  of  "precursor  front"  moving  ahead  of  the  actual  shock  front. 

The  width  of  this  front  can  be  conveniently  taken  as  that  range  of  £  over 
which  the  argument  of  the  erfc  function  changes  from  -1  to  +1.  From  the 
conditions  given  in  (lU)  and  (l5),  this  range  can  be  written  as 

M  «  i  -  i  -  4  <l6> 

It  is  seen  that  the  front  broadens  with  increasing  t,  and  for  any  t  is 
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sharper  for  smaller  D.  It  will  be  found  shortly,  however,  that  this  front 
is  so  broad  as  to  invalidate  the  idea  of  a  density  "front"  in  the  usual 
sense . 


In  order  to  obtain  a  more  precise  picture  of  this  behavior,  it  is 
convenient  to  introduce  dimensionless  parameters  into  Eq.  13.  If  we  write  D 
in  the  alternate  form,  D  =  v^L/3,  where  L  is  the  mean  free  path  for  electrons 
with  thermal  velocity  vq  in  the  neutral  gas  ahead  of  the  shock,  then  the 
arguments  of  the  functions  in  Eq.  13  may  be  put  into  the  form: 


IN  -3?v 


Ci- 

4ll>i 


=  AI 


with  the  dimensionless  parameters  T)  ,  rj  '  and  p  defined  by 


(i?) 


Tj  *i/ l  ,  y'.v.i./L  J  p.u-./o-.  (is) 


For  medium-weight  gases  (argon,  nitrogen,  etc.),  hypersonic  shocks  with  Mach 

numbers  of  the  order  of  10  have  velocities  of  about  U  x  10^  cm/sec,  and 

equilibrium  temperatures  of  the  order  of  10^  degrees  Kelvin.  Actually,  if 

we  assume  the  electron  temperature  to  be  equal  to  the  equilibrium  temperature 

of  the  gas,  then  for  hypersonic  shocks  in  ideal  gases,  v  and  v  are  each 

s  o 

approximately  proportional  to  the  Mach  number,  so  the  ratio  p  =  v  /v  is 

s  o 

—3 

independent  of  Mach  number  and  equal  to  about  6  x  10-J  .  On  the  other  hand, 
the  mean  free  path,  1,  can  differ  by  an  order  of  magnitude  from  one  gas  to 
another,  and  is  often  a  sensitive  function  of  electron  temperature.  With 
this  in  mind,  a  reasonably  informative  presentation  of  the  transient  behavior 
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described  by  Eq.  13  can  be  obtained  for  an  ideal  gas  by  setting  p  =  6  x  10"° 
and  plotting  log  (n/n^)  against  t)  for  several  values  of  the  normalized  shock 
front  position  r)  '  .  Such  a  family  of  curves  is  plotted  in  Fig,  k  „ 

The  curves  in  Fig.  1;  fail  to  suggest  the  existence  of  any  kind  of  well- 
defined  precursor  "front".  The  reason  for  this  may  be  seen  by  writing  the 
front  width  given  in  Eq.  16  in  terms  of  the  dimensionless  parameters  defined 
in  (18): 

=  A^/L  =  4  W/3L  *  4  ]7]73f'  (19) 

The  exponent  of  the  steady-state  decay  factor  e”^^  therefore  changes  over 
the  transition  region  by  an  amount 

3j>A9?  =  4  \if0p1 

This  means  that  the  width  of  the  transition  region,  even  for  very  small  values 
of  t ,  is  quite  a  bit  greater  than  a  decay  length  for  the  steady-state 
distribution,  so  the  effect  of  the  modulating  erfc  function  in  Eq.  13  is  spread 
over  a  large  part  of  the  precursor  distribution. 


HI®  Some  Implications  of  the  Transient  Behavior 

Although  the  precursor  density  does  not  possess  a  well-defined  front,  it 
is  possible  that  under  suitable  experimental  circumstances  the  behavior  shown 
in  Figo  U  could  give  rise  to  a  virtual  front  propagating  ahead  of  the  shock. 
The  existence  of  such  a  virtual  front  would  depend  on  the  nature  of  the 
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interaction  of  the  measuring  apparatus  with  the  precursor  density ,  and  could 
occur  under  several  experimental  conditions «  Let  us  list  the  characteristics 
of  some  probes  which,  if  placed  along  the  path  of  the  shock,  would  imply  the 
existence  of  a  precursor  front* 

1.  A  probe  fails  to  respond  until  the  electron  density  reaches  a 
certain  level * 

2.  A  probe  saturates  at  a  given  electron  density* 

3.  The  response  of  a  probe  changes  dis continuously  at  a  critical 
electron  density,, 

In  each  of  these  cases,  the  nature  of  the  response  of  a  probe  of  some  sort 
depends  critically  on  a  certain  electron  density,  and  it  is  the  propagation 
of  this  critical  level  that  constitutes  the  advancing  "front". 

A. )  Comparison  of  the  theory  with  some  experimental  results. 

The  critical  level  involved  in  case  1  might  be  simply  the  local  electron 
density  required  to  bring  the  response  of  the  probe  above  the  noise  level,  or 
some  other  threshold,  associated  with  the  detection  system*  This  could  have 
been  the  basis  of  the  "precursor  front"  observed  by  Weymann  using  induction 
probes. ^  He  claims  that  the  precursor  signal  had  a  pronounced  front,  but 
a  very  rapidly  growing  electron  density  rising  from  the  detection  limit  of 
his  apparatus  could  simulate  such  a  front.  Let  us  assume  that  this  was  the 
case,  and  attempt  to  match  the  propagation  of  a  level  of  constant  density  on 
Fig.  ii  with  the  x-t  diagram  of  a  measured  precursor  signal  front.  The  diagram 
we  will  use  is  Fig*  2  of  reference  1.  This  figure  has  been  reproduced  on 
Fig.  $  of  this  report. 
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The  conditions  under  which  the  measured  data  were  obtained  are  not 

provided,  but  they  may  be  guessed  from  information  given  elsewhere  in  the 

paper.  For  example,  in  Fig.  3  Weymann  gives  formative  time  lags. for  various 

diaphragms  bursting  into  either  3  or  6  mmHg  of  argon.  The  only  lags  as 

short  as  that  in  his  Fig.  2  occur  for  pressures  of  3  mmHg.  Talcing  this  as 

the  pressure  of  the  driven  gas,  we  find  in  his  Conclusions  estimates  of 

electron  temperature  (about  10^  degress  K),  and  diffusion  coefficient  (about 

1*  x  10  cm  /sec),  from  which  one  obtains  a  mean  free  path,  L,  of  about  0„2cm. 

Thus  we  can  establish  the  normalized  distances  tj  =  Z/L  and  i)'  =  v  t/L  in 

s 

our  Fig.  U  in  terms  of  the  measured  distances  in  his  Fig.  2.  This,  in  turn, 
makes  it  possible  to  look  for  a  level  of  constant  density  in  Fig.  U  whose 
propagation  is  consistent  with  the  signal  front  motion  Weymann  has  measured. 

We  choose  two  shock  front  positions  on  Fig.  f>,  say  (v  t),  =  h3  cm  and 

S  1. 

(v  t)  =  76  cm,  and  measure  the  distance  from  shock  to  signal  front  in  each 
s  •  *■ 

case.  We  find  that  the  signal  front  has  moved  relative  to  the  shock  front 
from  5^  =  62  cm  to  ^2  ~  cmJ  while  the  shock  front  itself  has  moved 
(76  -  li3)cm  =  33cm.  Using  L  =  0.2  cm  to  convert  to  our  normalized  coordinates, 
we  get 

V,  =  ijL  =  310  ,  y,  .  4»/l  * 

while 

-  A(uvO/(_  = 

With  these  data,  we  now  enter  Fig.  4  with  a  horizontal  line  whose  end¬ 
points  must  simultaneously  lie  on  the  vertical  lines  =  310,  t)  2  =  j?20  and 
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span  two  decay  curves  for  values  of  r)^  and  tj'  such  that  t}  '  -  tj'  =  165. 
The  line  drawn  on  Fig.  4  at  about  log(n/n^)  =  -11.5  approximately  satisfies 
this  requirement  by  lying  between  =  310  and  r\^  =  520  while  spanning  the 
decay  curves  for  =  30  and  rj'  =  200.  The  motion  of  a  level  of  constant 
density  along  this  line  is  plotted  in  Fig.  5  for  comparison  with  Weymann's 
measurements.  We  see  that  a  slightly  larger  formative  time  lag  is  obtained 
by  extrapolating  back  to  t  *  0,  and  the  first  measured  point  lies  below  the 
predicted  curve.  However,  the  model  assumes  the  instantaneous  creation  of 
and  ideal  shock  at  t=  0,  so  cannot  be  expected  to  hold  very  well  for  small 
values  of  t.  The  results  obtained  by  matching  at  the  same  points  for  dif¬ 
ferent  values  of  L  are  almost  identical  with  those  shown  in  Fig.  5  for 
1  =  0.2  cm.  A  larger  mean  free  path  would  allow  Weymann's  probes  to  have 
been  less  sensitive. 

It  is  of  interest  to  observe  that  as  t  -*  »  ,  the  velocity  of  the 
given  density  level  approaches  the  shock  velocity.  This  is  to  be  expected, 
the  transient  giving  way  to  the  fixed  steady-state  distribution  for  large 
times.  However,  over  the  range  for  which  the  experimental  data  are  given, 
we  may  assume  that  rj  >>  t)'  »  so  the  approximation  given  in  (15)  can  be  used 
for  the  precursor  density.  Clearly,  the  exponential  factor  dominates  the 
decay,  so  we  may  write 

°c  e.  (20) 

where  we  have  taken  K  +  v  t  -  (x  -  v  t)  +  v  t  =  x.  If  the  density  is  to 

o  SS 
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remain  constant,  so  must  the  exponent.  Thus  the  point  x  at  which  this 
constant  density  is  observed  must  vary  with  time  according  to 

x  -*c 


so  the  speed  of  this  level  in  the  laboratory  system  has  the  time-dependence 

(21) 


ctx  ^ 


37  -  « 


One  should  not  infer  from  this  equation  that  dx/dt  -*•  0  for  large  t.  For 
large  t,  the  approximation  t)  >>  t]  '  no  longer  holds,  and  we  must  find 
dx/dt  by  requiring  that  dn/dt  *  0  in  Eq.  13.  It  is  easily  shown  that  as 
t— ,  dn/dt  =  0  implies  dx/dt  -*  vg  as  shown  in  Fig.  Weymann  con¬ 
cluded,  presumably  from  several  sets  of  measurements,  that  the  velocity  of 
the  signal  front  in  shock  tube  coordinates  could  be  approximately  represented 
by  Eq.  21.  This  behavior  can  be  verified  for  the  motion  of  the  level  of 
constant  density  indicated  on  Fig.  it.  The  velocity  of  the  transient  "pre¬ 
cursor  front",  relative  to  the  shock  front,  may  be  written 

At  dfatVL  dl)'  >  (22> 

therefore  the  velocity  of  the  precursor  front  relative  to  the  laboratory 
system  becomes 


The  derivative  dr/dr) '  can  be  estimated  from  Fig.  k  for  various  values  of 
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•n'  (hence  t  ,  since  t  =  Lrj  '/v  )  along  the  line  of  constant  density 

s 

log(n/n  )  =  -11.6  when  L  =  0.2  cm.  A  plot  of  Eq.  23  is  given  in  Fig.  6, 

3? 

along  with  points  obtained  from  Weymann's  empirical  formula  (or  our  Eq.  21) 
by  normalizing  the  two  expressions  at  t  =  1E>0  sec.  The  agreement  over  the 
experimental  range,  t<15>0,  is  seen  to  be  good,  as  one  might  expect.  Our 
model  is  therefore  capable  of  an  a  priori  prediction  of  at  least  one  qualitative 
feature  of  the  experimental  results. 

The  arguments  presented  here  are  not  intended  to  "verify”  the  predictions 
of  the  simple  heuristic  theory  developed  in  this  report,  but  rather  to 
demonstrate  that  these  predictions  are  to  some  extent  consistent  with  the  only 
relevant  measurements  published  so  far.  However,  armed  with  a  bit  more  con¬ 
fidence,  we  will  go  on  to  discuss  the  implications  of  such  a  transient  pre¬ 
cursor  front  in  the  scattering  of  electromagnetic  waves. 

B.)  Effect  of  the  transient  on  the  scattering  of  electromagnetic  waves. 

If  a  plane  electromagnetic  wave  is  normally  incident  upon  the  shock, 
the  presence  of  the  precursor  electrons  ahead  of  the  front  can  be  expected 
to  affect  the  scattering  properties  of  the  shock.  Typically  in  problems 
of  this  sort,  there  are  many  parameters  to  be  considered.  Some,  like  the 
frequency  of  the  incident  wave  or  the  density  of  the  driven  gas,  may  vary 
over  wide  ranges.  Others,  like  the  shock  speed,  the  equilibrium  gas  tem¬ 
perature  and  the  electron  mean  free  paths,  are  not  precisely  known  and  must 
be  estimated.  Still  others,  like  the  model  parameters  n  and  T  ,  are  not 
known  at  all  and  must  be  guessed.  In  this  section  we  will  assume  driven  gas 
pressures  of  a  few  mmHg,  equilibrium  gas  temperatures  of  about  10^  degrees  K, 
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and  shock  speeds  of  about  h  x  10^  cm/sec.  The  precursor  electron  temperature 
will  be  taken  to  be  equal  to  the  temperature  of  the  hot  gas.  Electron  mean 
free  paths  ahead  of  the  shock  can  be  either  small,  as  in  air,  or  large,  as 
in  the  noble  gases,  the  ratio  being  ten  or  twenty  to  one.  A  reasonable  guess 
for  the  model  parameter  n^  is  very  difficult  to  justify,  and  therefore  requires 
more  extended  discussion. 

As  part  of  Weymann ' s  early  investigations,^  he  used  a  magnetic  probe  to 
obtain  an  estimate  of  the  precursor  density  far  ahead  of  the  shock.  He  con¬ 
cluded,  by  an  indirect  argument,  that  for  a  Mach  12  shock  in  limmHg  of  argon, 
the  precursor  density  is  about  10  /cur  at  a  distance  of  100  cm  ahead  of  the 
shock  front.  Later  measurements^  showed,  as  was  mentioned  in  the  Introduction, 
that  the  spatial  dependence  of  the  precursor  distribution  is  given  by  our  Eq.  1, 

which  implies  Eq.  2,  from  which  n  can  be  obtained  from  the  measured  density. 

6  2^  £ 

Using  D  -  3  x  10  cm  /sec  (for  ltfnmHg  argon)  and  v  =  I4.  x  1CT cm/sec,  we  find 

s 

12  3 

that  n^  :  6  x  10  /cm  .  This  seems  like  a  large  number.  If,  however,  the 
level  of  equilibrium  ionization  in  this  shock  is  calculated  from  the  information 
given  by  Resler,  et  al,  it  is  found  that  ng^  Z  2  x  101'Vcm^,  so  n^  is 
several  orders  of  magnitude  below  the  equilibrium  density  behind  the  onset 
point.  (It  is  also  of  interest  to  observe  that  in  order  to  make  the  electron 
energy-exchange  mechanism  work  in  a  low-voltage  arc  in  argon  at  a  few  mmHg, 
Druyvesteyn  found  that  electron  densities  of  the  order  of  10/cnr  were 
required. 

To  simplify  the  scattering  picture,  let  us  assume  that  the  decay  length, 

6  *  D/v  ,  of  the  steady-state  precursor  distribution  is  very  much  larger  than 

O 
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the  distance  from  the  shock  front  to  the  onset  point „  In  virtue  of  the  large 
electron  densities  behind  the  onset  point,  we  will  assume  this  region  to  be  a 
perfect  conductor,  and  move  it  up  to  the  shock  front  to  give  the  following  form 
to  the  scattering  obstacle: 

Perfect 

/l$u  ©c 


Pig.  7 

There  are  now  two  characteristic  numbers  which  determine  the  approx¬ 
imations  which  can  be  used  in  examining  the  scattering  properties  of  the 
shock.  One  of  these  is  the  maximum  plasma  frequency  encountered  in  the 
precursor  region;  this  is  cOpQ  =  (e  n^/me^)2  ~  2tt  x  10  rad/sec.  The  other 
is  a  characteristic  wavelength  equal  to  the  decay  length,  from  which  we  may 
define  a  frequency  cog  =  2tt  c/6  .  Unfortunately,  we  have  no  idea  of  how  n^ 
depends  on  the  type  or  pressure  of  the  driven  gas,  and  the  possible  range 


co  g  is  quite  large.  In  order  to  reduce  confusion,  let  us  assume  that 
cog<o0p0  ,  which,  if  our  various  hypotheses  are  correct,  will  be  the  case  for 
a  Mach  10-12  shock  in  a  few  mmHg  of  argon.  In  this  case  the  division  of  the 


frequency  spectrum  will  .look  something  like 


Figo  8 
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It  should  be  kept  in  mind  that  we  are  not  dealing  with  a  stationary 
scattering  problem*  The  shock  is  in  motion*  and  the  reflected  wave  carries 
knowledge  of  this  motion  through  a  time- varying  phase  factor*  The  scattering 
problem  should  properly  take  into  account  both  the  time  and  space  variations 
of  the  electron  distribution!  an  analysis  in  terms  of  the  continuous  frequency 
spectrum  of  the  scattered  wave  would  be  appropriate*  However,  the  shaky  struc¬ 
ture  of  hypothesis  and  inference  erected  thus  far  can  hardly  justify  such  a 
detailed  analysis*  It  will  be  adequate  for  our  purpose  to  consider  the  shock 
front  fixed  when  calculating  the  effect  of  the  precursor  electrons  on  the  in¬ 
cident  and  reflected  waves*  The  "ether  drag"  associated  with  the  motion  of 
the  precursor  will  give  perturbations  of  opposite  sign  for  the  two  waves,  and 

p 

can  therefore  be  expected  to  contribute  corrections  of  order  (v  /c)  *  Shock 

s 

speeds  being  what  they  are,  this  is  negligible*  Accordingly,  any  change  in 

the  Doppler  effect  will  be  considered  to  arise  from  changes  of  the  precursor 
distribution  relative  to  the  moving  shock,  rather  than  to  a  fixed  observer. 

Scattering  at  various  frequencies  below,  between  and  above  the  reference 
points  in  Fig.  8  will  now  be  discussed  in  a  qualitative  way* 

1 . )  Low-Frequencies :  to  «  . 

At  these  low  frequencies,  the  steady-state  precursor  distribution  extends 
only  a  fraction  of  an  incident  wavelength  from  the  shock  front,  and  may  be  re¬ 
placed  by  an  equivalent  plane  sheet  of  thickness  6  and  electron  density  n^. 
Whatever  its  properties,  such  an  electrically  thin  layer  cannot  be  expected  to 
change  the  reflecting  properties  of  the  region  over  which  it  is  placed  in  any 


AF-4561/ll 


23 


significant  way.  This  can  easily  be  verified.  For  the  transient  case*  the 
effect  of  a  perturbation  that  grows  in  a  reasonable  time  interval  from  zero 
to  an  insignificant  value  is  still  insignificant,  We  will  conclude,  then, 
that  at  sufficiently  low  frequencies,  the  electron  precursor  will  produce  no 
observable  effects. 


2.)  High-Frequencies ;  co  »  ®  . 

For  sufficiently  high  frequencies,  the  scattering  amplitude  can  be  cal¬ 
culated  in  the  geometrical  optics  approximation.  Consider,  for  the  moment, 
a  stationary  scattering  problem.  Assuming  a  time-dependence  exp  (-j  cot) ,  the 
electric  field  in  the  presence  of  a  space-varying  electron  density  n(x)  must 


satisfy  the  equation 


|  +  (C  -p?c«>)E.  =  O 


where 


C  =  ,  &«)  *  e* ncx)/g0*vca(i  -vjui)  . 


2  2 

Under  the  conditions  given  above,  8  <<  k  and  (dB  /B  dx)  «  k  ,  so  the  con- 

r  ^  F  r  W 

ditions  for  the  application  of  the  geometrical  optics  approximation  are  satis¬ 
fied,  and  we  may  write 


fc.  (x)  OC  £ 


+ cSx 


To  avoid  needless  complication,  it  may  be  reasonably  assumed  that  v  /co  «  1, 

v 

2 

making  p  a  real  positive  number 0 

P  '  -  -  ---  
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It  will  be  helpful  to  consider  first  the  degenerate  case  in  which  there 
are  no  precursor  electrons.  A  reference  plane  (which  could  represent  a  re¬ 
ceiving  antenna)  will  be  established  at  xq}  and  the  shock  front  will  be  loca¬ 
ted  at  x  <  x  «  as  in  the  figure  belows 
s  o* 


-<i0 e.x+w-6) 
R<fei -e*ee  ^ 

I  |  I  I 


X-O 


Xs 


Xo 


Pig.  9 

The  incident  wave  at  xQ  is  simply  Einc  =  exp(-j(koxQ+wfc)  ).  But  this  wave 
moves  past  xq  up  to  xgJ  where  it  is  reflected  (phase  increase  of  180°)  and 

p 

sent  back  toward  xq.  Using  Eq„  26  with  =  0,  we  can  calculate  the  phase 

accumulation  factor  in  the  round  trip  from  x  to  x  to  x  0  Adding  jtr  to 

o  s  o 

account  for  the  reflection  by  a  perfect  conductor 5  and  multiplying  by  the 

incident  wave  at  x  ,  we  get  the  reflected  wave  at  x  in  the  form 

o 

£  <Sx  + 

E.r^,60  =  -  e  **  a  ,  (27) 


which  can  be  rewritten  as 

■^(^oXo  —  <js£  — 

E^|(Xo)  a  -a  .  (28) 
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Let  us  now  reintroduce  the  motion  of  the  shock  by  setting  x  <*  v_t,  and 

s  s 

define  the  generalized  frequency  of  a  wave  at  a  fixed  point  as  the  time 
rate  of  change  of  its  phase 3  that  is, 


/ 

LA 


ct 

^  It  ‘ 


(29) 


For  the  reflected  wave  given  in  Eq.  28,  the  phase  is 

$  *  “iCWXo  -wt  - 
so  its  frequency  is 


t 

wre4l 


(A  + 


CO  +  ~ 

/ 


(30) 


thus  displaying  the  Loppler  shift  produced  by  the  motion  of  the  shock. 

The  perturbation  introduced  by  the  steady-state  precursor  distribution 
can  be  calculated  by  writing  the  exponent  in  Eq.  26  in  the  approximate  form 


4^.(1“^)/^)^ 


(31) 


The  first  term  contributes  the  phase  component  we  have  already  discussed.  In 

view  of  our  earlier  discussion,  the  second  tern  is  to  be  calculated  relative 

to  the  moving  shock.  Therefore,  transforming  to  the  shock -based  coordinates, 

o 

and  assuming  that  xq  is  sufficiently  far  away  that  {3^  (xq)  «  0  for  all  times 
of  Interest,  we  may  write  the  precursor  contribution  to  the  phase  $  as. 
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This  is  simply  a  fixed  phase  shift  in  the  reflected  wave,  and  could  therefore 
contribute  nothing  of  observable  interest* 

2 

Under  transient  conditions,  however,  P  is  growing  from  effectively  zero 
to  its  steady-state  value*  Since  the  effective  dielectric  constant  of  the 
electron  distribution  is  less  than  unity,  the  growth  of  the  precursor  will  be 
accompanied  by  a  decrease  with  time  of  the  total  phase*  This  will  result  in 
an  increase  in  the  Doppler  shift  inferred  from  the  reflected  wave.  The  time 
dependence  of  that  part  of  the  phase  due  to  the  precursor  can  be  introduced 
by  altering  the  first  equality  of  (32)  to  read 


where  the  integration  is  again  performed  in  shock  coordinates, 
of  the  reflected  wave  now  becomes 

°o 


The  frequency 


UY  ,,  =  a  =  w  +  ^V.vr*  + 

J  d-t 


I 


9*. 


(3U) 


2 

The  time  derivative  of  p  may  be  found  by  using  n(£,t)  from  Eq*  13*  The  time 
differentiation  and  subsequent  ^-integration  are  straightforward,  though  com¬ 
plicated,  and  we  will  record  only  the  results 


0$) 
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For  an  upper  limit  of  the  effect  produced  in  this  way,  let  us  calculate  A©' 

P 

for  k  =  P  ,  v  «*  ipcl.O'’ em/sec  and  D  «  IpdO^cm^/sec,  (the  driven  gas  is 
o  po?  s 

.  t 

argon  at  about  3mmHg.)  The  unperturbed  Doppler  shift  is  Ago  =  2k  v_*  Below 

IX  c  s 


We  see  that  the  perturbation  here  is  at  most  a  small  fraction  of  the  normal 
Doppler  effect o 


3.) 


Intermediate-Frequencies:  ©_<©<©  • 

. o  po 

Over  a  certain  restricted  range  of  frequencies,  it  is  possible  that  the 


transient  behavior  of  the  precursor  may  produce  a  more  pronounced  effect  than 


in  the  two  extreme  cases  just  discussed,,  For  ©  in  the  range  given  above,  one 


can  still  assume  that  the  precursor  extends  over  several  incident  wavelengths 
(©g  <  ©),  but  for  some  value  of  K  >  0,  the  coefficient  of  E  in  Eq„  2k  vanishes. 
This  means  that  the  effective  dielectric  constant  of  the  precursor  vanishes  at 


some  point  ahead  of  the  shock,  and  the  incident  wave  is  perfectly  reflected  at 
that  point.  For  a  given  frequency  <*s  this  point  is  defined  by  the  equation 
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AF- 


or,  solving  for  n  in  shock  coordinates  (since  x„  =  x  +  S_), 

r  s  r 

n(ir)  -  ^  cJ'/z*'  ■ 


(36) 


In  the  steady-state,  the  incident  wave  is  simply  reflected  from  a  point  at 
the  fixed  distance  ahead  of  the  actual  shock  front,  and  the  phase  of  the 
reflected  wave  is  changed  by  a  fixed  amount,. 

During  the  transient  period  of  precursor  accumulation,  however,  the 
point  of  reflection,  moves  forward  from  the  actual  shock  fronto  That  is, 
since  this  point  corresponds  to  the  critical  density  defined  by  the  right  side 
of  Eq.  36,  the  propagation  of  this  density  level,  as  discussed  at  the  beginning 
of  the  section,  will  produce  a  corresponding  movement  of  the  point  of  reflec¬ 
tion.  The  appropriate  density  can  be  represented  by  a  horizontal  line  on  Fig.  I4., 
and  the  velocity  of  the  reflection  point  will  resemble  the  precursor  "front" 
velocity  shown  in  Fig.  6.  This  will  show  up  as  an  enhanced  Doppler  shift  at 
the  beginning  of  the  shock. 

We  can  make  these  observations  a  bit  more  quantitative  by  using  the  WKB 
approximation^^- ^  in  solving  Eq„  2lu  We  begin,  as  in  the  high-frequency  case, 
by  considering  a  stationary  scattering  problem.  Taking  xq  in  Fig.  9  suffici¬ 
ently  far  away  from  the  shock,  the  asymptotic  form  of  the  local  turning-point 
solution  may  be  useds^^ 


r 

Le 


-i(+r-rK)~ 

e 


(37) 


where 
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NAT  = 


-  Op  (<) 


The  constant  A  is  to  be  determined  by  the  reqairement  that  the  electric  field 

E  at  x  be  the  sum  of  the  incident  wave  fend  a  reflected  wave.  One  way  to  find 
o 

such  a  constant  is  to  add  and  subtract  kQ  in  the  integrand  of  Eq.  38,  writing 
w  in  the  form 


vr  »  V  =  k8(xe-Xr)  - 


Where 


4fx)  =  j~l?o  -  4  . 


Reintroducing  the  time  dependence  exp(-j©t)  into  Eq.  37,  we  obtain. 


E  C/oj-6)  —  E  (yj  6 


--j^oXe,  +^oXr 


+  e. 


The  incident  wave  at  xq  must  be  exp(-j  (kQxo  +  <ot)  ),  and  can  be  identified  in 
the  second  term  in  the  brackets  above.  We  choose  A  in  such  a  way  that  the 
incident  wave  is  left  standing  alone,  by  taking 


.  ^E5  e~^lu'p  ~ -5  £^fy>Ax 
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This  accomplishes  the  required  decomposition  of  E(xQ,t); 

E  (*»£)  -  e  +  ^  je  .  0»3) 

The  reflection  coefficient  is  the  expression  in  brackets.  Having  isolated  the 
reflected  wave,  the  various  terms  in  its  phase  will  now  be  recombined,  giving 

§teHM  =  ~j  (i.*.  +  <»t  Oil.) 

^  r 

The  integral  in  this  equation  extends  from  the  reflection  (or  turning)  point 

2 

x^  to  the  far-removed  observation  point  xq.  Accordingly,  is  different  from 

zero  over  only  a  small  part  of  the  range  of  integration,  and  a  reasonable  ap- 

2 

proximation  can  be  expected  if  we  simply  ignore  B  and  write 

Jr 

f*’ 

\ J  a?) 

Equation  UU  then  becomes 


^  ^1  (-6)  «  — ; j  jj-  W0Xo  4-  Got  +  3ke,  4  j  (U6) 

where  we  have  written  x,,  »  x_  +  £  (t)  =  v  t  +  ?_(t)«  According  to  (29),  the 

r  s  r  sr 

frequency  of  the  reflected  wave  becomes 


( 


w  +  a.  k0\r5  +  ■at,, 


jj/*) 

It  \ 


Off) 


AF-it.561/11 


31 


The  time  dependence  of  ?r(t)  is  most  easily  obtained  by  reference  to  Fig.  It. 

A  given  frequency  defines  an  electron  density  through  Eq.  36,  and  thus  a 
horizontal  line  on  Fig.  U.  The  velocity  of  propagation  of  this  level  relative 
to  the  shock  front  can  then  be  found  numerically  in  terns  of  increments  dr)  and 
dr]  .  (See  Eq.  22  and  the  discussion  following  it.) 

If  the  electron  density  upstream  of  the  reflection  point  is  taken  into 
account  in  (I|5),  the  effect  is  to  produce  a  slight  additional  time  varying  de¬ 
crease  in  the  phase.  This  will  augment  the  Doppler  effect  still  further,  since 
the  smaller  the  phase,  the  closer  the  reflection  point  appears  to  be  to  the 
observation  point  xq.  We  will  not  bother  with  this  small  effect,  other  than  to 
observe  that  as  a  result,  the  approximation  (U7)  is  a  lower  limit  to  the  en¬ 
hanced  Doppler  shift o 

The  magnitude  of  the  shift  will  depend  upon  the  size  of  to  compared  to 
The  value  of  n^  inferred  from  Weymann's  measurements  is  np  ~  6xl012/cm^, 
so  let  us  be  as  optimistic  as  possible  and  assume  that  such  a  value  is  not  un¬ 


realistic.  If  we  take  a  few  representative  frequencies,  2  kmc,  5  kmc  and  8  kmc, 

we  find  the  critical  densities  defined  by  Eq.  36  are  8.3x10“^  n  ,  £.2xlCf2  n 

P  P 

and  0.13  n^,  respectively.  These  numbers  define  the  dotted  horizontal  lines 
shown  on  Fig.  lu  The  effective  shock  velocity,  inferred  from  the  Doppler  shift 
of  the  reflected  wave,  can  be  estimated  from  Eq.  23  for  each  of  these  frequen¬ 


cies. 


The  results  are  plotted  in  Fig.  10  as  v  vs.  t,  where  the  driven  gas 


is  again  3  mmHg  argon.  We  see  that  the  lower  the  frequency,  the  greater  the 


effect.  It  should  be  remembered,  however,  that  the  WKB  approximation  requires 


that  the  electron  density  should  not  decay  too  rapidly  over  an  incident  wave¬ 


length.  In  fact,  this  is  already  violated  for  the  lower  frequencies  in  the 
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example  above. 

It  is  unlikely  that  any  of  the  effects  discussed  here  could  be  observed 
in  air,  the  main  reason  being  that  the  ratio  vg/D  is  so  much  larger  in  air 
than  in  an  ideal  gas  like  argon.  The  equilibrium  temperature  of  the  shock 
heated  gas  will  be  less  for  air  than  for  an  ideal  gas  at  the  same  shock  speed. 
More  important,  the  collision  probability  for  electrons  at  the  usual  shock 
temperatures  is  more  than  an  order  of  magnitude  greater  in  air.  This  means 
that  the  mean  free  path  is  so  small  that  any  transient  effects  that  might 
occur,  can  have  a  duration  of  only  a  few  microseconds.  For  example,  in 
Fig.  6,  the  reduction  of  mean  free  path  by  a  factor  l/lO  will  divide  the 
t-scale  by  10  and  everything  would  be  over  after  about  l£  p  sec.  Moreover, 
the  approximation  on  which  Fig.  10  is  based  would  be  completely  invalid, 
since  D/v  would  now  be  only  about  1  cm,  while  the  shortest  incident  wave- 
length  is  almost  U  cm. 


IV.  Summary  and  Conclusions 

How  much  confidence  can  be  placed  in  the  electron  diffusion  hypothesis 
underlying  the  results  obtained  here?  The  answer  to  this  question  is  related, 
in  part,  to  the  results  themselves.  Actual  electrons,  as  well  as  electrical 
effects  allegedly  due  to  them,  have  been  detected  far  ahead  of  shock  fronts, 
and  the  most  compelling  explanation  for  these  observations  is  that  the  electrons 
got  there  by  diffusion.  This  hypothesis  is  immediately  called  into  question  by 
a  vague  expectation  that  strong  space-charge  fields  would  prohibit  any  appre- 
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ciable  charge  separations  But  without  an  adequate  theory  of  this  admittedly 
complex  physical  situation,  there  is  no  more  reason  to  reject  than  to  accept 
the  diffusion  hypothesis.  We  have  attempted  to  give  credence  to  the  idea  of 
electron  diffusion  through  the  shock  front  by  suggesting  a  few  mechanisms 
that  are  thought  to  operate  in  different,  but  somehow  analogous,  situations. 

In  the  end,  we  are  led  to  accept  the  hypothesis,  at  least  provisionally,  and 
to  ask  what  its  Implications  might  be,  and  whether  these  implications  are  con¬ 
sistent  with  the  existing  experimental  data.  The  model  we  have  chosen  to 

study  is  heuristic,  in  the  sense  that  it  emphasizes  the  particular  features  of 
the  hypothesis  we  wish  to  examine,  and  leads  to  fairly  clear  predictions  of 
behavior  which  can  be  measured.  Comparing  these  predictions  with  the  data 
already  available,  we  discover  the  following : 

1. )  The  spatial  dependence  of  the  predicted  steady-state  electron  dis¬ 

tribution  agrees  with  both  the  elementary  theory  and  the  measured 
results. 

2. )  The  transient  behavior  of  the  precursor  distribution  obtained  from 

the  model  suggests  the  possibility  of  observing  a  "precursor  front" 
under  certain  conditions.  Allowing  a  reasonable  interpretation  of 
the  way  in  which  the  observed  precursor  front  was  measured,  it  is 
found  that  the  predictions  of  the  model  are  consistent  with  the  ex¬ 
perimental  results.  In  fact,  the  empirical  time-dependence  of  the 
observed  precursor  front  velocity  is  given  exactly. 

Furthermore, 

3. )  It  is  suggested  that  the  transient  should  produce  a  relatively  pro- 
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nounced  effect  on  the  scattering  of  electromagnetic  waves  in  a 
certain  frequency  range,  (for  suitable  gases  and  pressures.) 

An  experiment  is  being  planned  to  check  this  behavior. 

In  short,  these  results  have  added  a  few  more  pieces  to  the  pattern  of 
plausible  inference  which  presently  supports  the  electron  diffusion  hypothesis. 
If  the  scattering  experiment  mentioned  above  should  prove  successful,  further 
supporting  evidence  would  be  provided. 

Theory  and  experiment  agree  that  the  dramatic  precursor  effects  are  to  be 
found  in  the  noble  gases.  For  shocks  in  air,  the  lower  gas  temperatures  and 
much  smaller  mean  free  paths  combine  to  limit  sharply  both  the  spatial  range 
of  the  precursor  distribution  and  the  duration  of  the  transient.  However, 
experiments  using  the  noble  gases  are  likely  to  provide  additional  information 
about  the  electrical  characteristics  of  the  shock  front  region. 
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FIG.  4  BEHAVIOR  OF  PRECURSOR  ELECTRON  DENSITY  AS  A  FUNCTION 
OF  £  AND*,  (p  =  Vs  /V0  =6X10'*) 


009 


FIG.  5  COMPARISON  OF  MEASURED  SIGNAL  FRONT  WITH  THAT  OBTAINED  FROM  FIG.  4 
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FIG.  10  EFFECTIVE  SHOCK  VELOCITY  FROM  DOPPLER  SHIFT 
FOR  FREQUENCIES  2  kmc,  5 kmc,  8  kmc. 
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Library  2 

National  Bureau  of  Standards 
Boulder  Laboratories 
Boulder,  Colorado 

National  Bureau  of  Standards  1 

U.  So  Department  of  Commerce 
Washington  23.  D.  C. 

Attn;  Mr.  A.  G.  McNish 

National  Bureau  of  Standards  1 

U.  3.  Department  of  Commerce 
Washington  23*  D.C. 

Attn;  Gustave  Shapiro,  Chief 

Engineering  Electronics  Section 
Electricity  and  Electronics  Division 

Director  7 

National  Security  Agency 
Washington  23,  D„C. 

Attns  R/D  (331) 

Office  of  Technical  Services  2 

Department,  of  Commerce 

Washington  23,  D„C0 

Attn;  Technical  Reports  Section 

AFCRL,  Office  of  Aerospace  Research  (CRRELT)  10 
La  G.  Hanscom  Field 
Bedford,  Massachusetts 

AFCRL,  Office  of  Aerospace  Research  (CRRD)  3 

Attn:  Carlyle  J.  Sletten 
L.  G.  Hanscom  Field 
Bedford,  Massachusetts 

AFCRL,  Office  of  Aerospace  Research  (CRRD)  2 

Attn:  Contract  Tiles 
Lo  G.  Hanscom  Field 
Bedford,  Massachusetts 
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Code 

M  $h 

M  5? 

G  116 

N  1 

N  16 

N  27 

N  29 

N  30 

N  37 


Organization 


Hq„  AFCCDD  (CCRRC,  Capt ,  John  J.  Hobson) 
L.  G.  Hans com  Field 
Bedford,  Massachusetts 

AFCCDD 

L.  G„  Hanscom  Field 
Bedford,  Massachusetts 

NASA  (Dr.  R.  V.  Hess) 

NASA  Langley  Research  Center 
Hampton,  Virginia 

Director,  Avionics  Division  (AV) 

Bureau  of  Aeronautics 
Department  of  the  Navy 
Washington  25,  D.C. 

Commander 

U.  S.  Naval  Air  Missile  Test  Center 
Point  Mugu,  California  . 

Attn:  Code  366 

Librarian 

II.  S.  Naval  Postgraduate  School 
Monterey,  California 

Director 

U.  S.  Naval  Research  Laboratory 
Washington  25,  D.C. 

Attn:  Code  2027 

Dr.  J.  I.  Bohnert,  Code  5210 
U.  S.  Naval  Research  Laboratory 
Washington  25,  D.C. 

Chief  of  Naval  Research 
Department  of  the  Navy 
Washington  25,  D.C. 

Attn:  Code  h27 


No.  of  Copies 

1 

1 

1 

1 

1 

1 

2 

1 

1 
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Code 

N  48 

N  85 

N  86 

N  88 

N  91 

N  lUl 

N  73 

I  1 

I  388 


Organization 


Commanding  Officer 
U.  S„  Haval  Air  Development  Center 
Johnsville,  Pennsylvania 
Attn:  NADC  Library 

Commanding  Officer  and  Director 

U„  S„  Navy  Electronics  Laboratory  (Library) 

San  Diego  52,  California 

Chief,  Bureau  of  Ordnance 
Department  of  the  Navy 
Surface  Guided  Missile  Branch 
Washington  25,  D.C. 

Attn:  Code  ReSl-e 

Department  of  the  Navy 
Bureau  of  Aeronautics 
Technical  Data  Division,  Code  i*106 
Washington  .25,  D.C0 

Commander 

U„  S„  Naval  Air  Test  Center 
Patuxent  River,  Maryland 
Attn:  ET-315,  Antenna  Branch 

ARDC  Regional  Office 
c/o  Department  of  the  Navy 
Room  4549,  Munitions  Building 
Washington  25,  D.C. 

Office  of  Naval  Research 
Branch  Office,  London 
Navy  100,  Box  39 
Fo  Po  0.  New  York,  New  York 

Airborne  Instruments  Laboratory,  Inc. 
Division  of  Gutler  Hammer 
Walt  Whitman  Road 
Melville,  L,  I.  New  York 
Attn:  Library 

Air com  Inc, 

351+  Main  Street 
Winthrop,  Massachusetts 


of  Copies 

1 

1 

1 

1 

1 

1 

10 

1 

1 


AF-4561/H 


List  S-S  -  Page  7 


Code  Organization 


I  37h  ACF  Industries*  Inc. 

Southeast  Comer*  5 2nd  Avenue  and 
Jackson  Street 

Bladensburg*  Prince  George’s  County 

Maryland 

Attn;  Librarian 

I  20 5  Battelle  Memorial  Institute 

5o5  King  Avenue 
Columbus  1*  Ohio 

Attn;  Wayne  E„  Rife*  Project  Leader 

Electrical  Engineering  Division 

I  8  Bell  Aircraft  Corporation 

Post  Office  Box  One 
Buffalo  5*  New  York 
Attn;  Eunice  P„  Hazelton*  Librarian 

I  2+69  Bell  Telephone  Laboratories 

Murray  Hill 
New  Jersey 

I  13  Bell  Telephone  Laboratories*  Inc. 

Whippany  Laboratory 
Whippany*  New  Jersey 
Attn;  Technical  Information  Library 

I  9l;0  Aerospace  Corp„ 

Box  95085 

Los  Angeles  !£*  California 
Attn;  Library 

I  2I4.7  Bendix  Radio  Division 

Bendix  Aviation  Corporation 
E.  Joppa  Road 
Towson  U*  Maryland 
Attn;  Dr0  D0  M,  Allison*  Jr. 

Director  Engineering  &  Research 

I  2I4.8  Bjorksten  Research  Laboratories,  Inc. 

P.  0.  Box  265 
Madison*  Wisconsin 
Attn;  Librarian 


No,  of  Copies 

1 

1 

1 

1 

1 

1 

1 
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Code 

Organisation 

No.  of  Copies 

I  2h9 

Boeing  Airplane  Company 

Pilotless  Aircraft  Division 

P.  0.  Box  3707 

Seattle  2li,  Washington 

Attn;  R.  R«  Barber,  Library  Supervisor 

2 

I  892 

Boeing  Airplane  Company 

7755  East  Marginal  Way 

Seattle  2hs  Washington 

Attn;  Mr.  D.  L.  Scidmore 

1 

I  253 

Chance  Fought  Aircraft,  Inc. 

93lU  West  Jefferson  Street 

Dallas,  Texas 

Attn;  A.  D.  Pattuilo,  Librarian 

1 

I  1*70 

Chu  Associates 

P.  0.  Box  38? 

Whitcomb  Avenue 

Littleton,  Massachusetts 

1 

I  918 

Collins  Radio  Company 

855  35th  Street,  N.E. 

Cedar  Rapids,  Iowa 

Attn;  Dr.  R.  L„  McCreary 

1 

I  919 

Collins  Radio  Company 

1200  North  Alma  Road 

Richardson,  Texas 

Attn;  C.  D.  Tipton 

1 

I  126 

Convair,  A  Division  of  General  Dynamics  Corp.  1 

Fort  Worth,  Texas 

Attn;  K,  G„  Brown. 

Division  Research  Librarian 

I  25U 

Convair,  A  Division  of  General  Dynamics 
2165  Pacific  Highway 

San  Diego  12,  California 

Attn;  Mrs.  Dora  B.  Burke 

Engineering  Librarian 

Corp.  1 
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Code 


Organization 


No.  of  Copies 


I  25 

i  255 

I  28 

I  257 

I  258 

I  259 

I  187 

I  ia5 


Cornell  Aeronautical  Laboratory,  Inc,  1 

4455  Genesee  Street 
Buffalo  21,  New  York 
Attn:  Librarian 

Dalmo  Victor  Company  1 

A  Division  of  Textron,  Inc. 

1515  Industrial  Way 
Belmont,  California 

Attn:  Mary  Ellen  Addems,  Technical  Librarian 

Dome  and  Margolin,  Inc.  1 

29  New  York  Avenue 

Westbury,  Long  Island,  New  York 

Douglas  Aircraft  Co.,  Inc.  1 

827  Lapham  Street 
El  Segundo,  California 
Attn:  Engineering  Library 

Douglas  Aircraft  Company,  Inc.  1 

3000  Ocean  Park  Boulevard 
Santa  Monica,  California 
Attn:  Peter  Duyan,  Jr. 

Chief,  Electrical/Electronics  Section 

Douglas  Aircraft  Company,  Inc.  1 

2000  North  Memorial  Drive 
Tulsa,  Oklahoma 

Attn:  Engineering  Librarian,  D-250 

Electromagnetic  Research  Corporation  1 

5001  College  Avenue 
College  Park,  Maryland 
Attn:  Mr.  Martin  Katzin 

Electronica  Communication  1 

1830  York  Road 
Timonium,  Maryland 
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Code 

I  2 99 

I  262 

I  11*7 

I  26U 

I  266 

I  269 

I  793 

I  893 


Organization 


No»  of  Copies 


Electronic  Specialty  Company  1 

E>121  San  Fernando  Road 
Los  Angeles  39,  Cal±f©xma 
Attn;  Donald  L.  Margerum 

Chief  Engineer,  Radiating  Systems  Division 

Emerson  Electric  Mfg.  Co.  1 

8100  West  Florissant  Avenue 

St.  Louis  21,  Missouri 

Attn;  Mr.  E„  R.  Breslin,  Librarian 

Emerson  Radio~Phonograph  Corp.  1 

Emerson  Research  Laboratories 
701  Lamont  Street,  N.W. 

Washington  10,  D.C. 

Attn:  Mrs.  R„  Corbin,  Librarian 

Fairchild  Aircraft-Missiles  Division  1 

Fairchild  Eng.  and  Airplane  Corp. 

Hagerstown,  Maryland 
Attn;  Library 

Federal  Telecommunication  Laboratories  1 

500  Washington  Avenue 
Nutley  10,  New  Jersey 
Attn;  Technical  Library 

General  Electric  Company  1 

Electronics  Park 
Syrcause,  New  York 
Attn;  Documents  Library 

B„  Fletcher,  Building  3-ll*3A 

General  Electric  Company  1 

Missile  and  Space  Vehicle  Department 

3198  Chestnut  Street 

Philadelphia ,  Pennsylvania 

Attn:  Documents  Library 

General  Electric  Company  1 

37^0  D  Street 

Philadelphia  2l*,  Pennsylvania 
Attn:  Mr.  H.  G„  Lew 

Missile  and  Space  Vehicle  Department 
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Code 

I  270 

I  ii8 

I  UU8 

I  737- 

I  207 

I  % 

I  265 

I  230 


Organization 


General  Precision  Laboratory,  Inc . 

63  Bedford  Road 
Pleasantville,  New  York 
Attn :  Lib  rarian 

Goodyear  Aircraft  Corp. 

1210  Massillon  Road 

Akron  1 5,  Ohio 

Attn:  Library,  Plant  G 

Granger  Associates 

Electronic  Systems 

97k  Commercial  Street 

Palo  Alto,  California 

Attn:  John  V.  N.  Granger,  President 

The  Hallicrafters  Co. 

5th  and  Kostner  Avenues 
Chicago  2k t  Illinois 
Attn:  Henri  Hodara,  Head  of 
Space  Communication 

Hughes  Aircraft  Company 
Antenna  Department 
Building  12,  Mail  Station  271k 
Culver  City,  California 
Attn:  Dr,  W.  H.  Kummer 

Hughes  Aircraft  Company 
Florence  Ave.  and  Teale  Streets 
Culver  City,  California 
Attn:  Louis  L,  Bailin 

Manager,  Antenna  Department 

ITT  Laboratories 
3700  East  Pontiac  Street 
Fort  Wayne  1,  Indiana 
Attn:  Technical  Library 

Jansky  and  Bailey,  Inc. 

1339  Wisconsin  Avenue,  N.W. 
Washington  7,  D.C. 

Attn:  Mr,  Delmer  C.  Ports 


No.  of  Copies 

1 


1 


1 


1 


1 


1 


1 


1 
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Code 


Organization 


No,  of  Copies 


I  21 +1 


I  279 


I  1+68 


I  136 


I  280 


I  282 


I  116 


I  61+8 


Dr.  Henry  Jasik,  Consulting  Engineer  1 

298  Shames  Drive 

Brush  Hollow  Industrial  Park 

Westbury,  New  York 

Lockheed  Aircraft  Corporation  1 

2555  N.  Hollywood  Way 

California  Division  Engineering  Library 

Department  72-25,  Plant  A-l,  Building  63-1 

Burbank,  California 

Attn:  N.  C.  Harnois 

Lockheed  Aircraft  Corporation  1 

Missile  Systems  Division  Research  Library 

Box  50l+,  Sunnyvale,  California 

Attn:  Miss  Eva  Lou  Robertson,  Chief  Librarian 

The  Martin  Company  1 

P.  0.  Box  179 

Denver  1,  Colorado 

Attn:  Mr.  Jack  McCormick 

The  Martin  Company  1 

Baltimore  3,  Maryland 
Attn:  Engineering  Library 
Antenna  Design  Group 

McDonnell  Aircraft  Corporation,  Dept.  661+  1 

Box  5l6 

St.  Louis  66,  Missouri 
Attn:  C.  E,  Zoller 

Engineering  Library 

Melpar,  Inc.  1 

3000  Arlington  Boulevard 

Falls  Church,  Virginia 

Attn:  Engineering  Technical  Library 

The  Mitre  Corporation  1 

21+1+  Wood  Street 

Lexington  73,  Massachusetts 

Attn:  Mrs.  Jean  E.  Claflin,  Librarian 
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Code 


Organization 


No.  of  Copies 


I  93k 


I  61a 


I  281+ 


I  285 


I  82 


I  286 


I  287 


I  225 


Motorola,  Inc .  1 

Phoenix  Research  Laboratory 

3102  N.  56th  Street 

Phoenix,  Arizona 

Attn:  Dr.  A.  L.  Aden 

National  Research  Council  1 

Radio  &  Electrical  Engineering  Division 
Ottawa,  Ontario,  Canada 
Attn:  Dr.  G.  A.  Miller,  Head 
Microwave  Section 

North  American  Aviation,  Inc.  1 

12211+  Lakewood  Boulevard 
Downey,  California 

Attn:  Technical  Information  Center  ( U95— 12 ) 
Space  and  Information  Systems  Division 

North  American  Aviation,  Inc.  1 

Los  Angeles  International  Airport 

Los  Angeles  1+5,  California 

Attn:  Engineering  Technical  File 

Northrop  Corporation  1 

Norair  Division 
1001  E.  Broadway 
Hawthorne,  California 

Attn:  Mr.  E.  A,  Freitas,  Library  Dept.  311+5 

Page  Communications  Engineers,  Inc.  1 

2001  Wisconsin  Avenue,  N.W. 

Washington  7,  D„C. 

Attn:  (Mrs.)  Ruth  Temple,  Librarian 

Philco  Corporation  1 

C  and  Tioga  Streets 
Philadelphia  3l+,  Pennsylvania 
Attn:  Mrs.  Dorothy  B.  Collins 
Research  Librarian 

Pickard  and  Burns,  Inc.  1 

21+0  Highland  Avenue 
Needham  9k,  Massachusetts 
Attn:  Dr.  Richard  H.  Woodward 
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Code 


Organization 


No.  of  Copies 


I  288 


I  232 


I  289 


I  91U 


I  290 


I  291 


I  757 


I  789 


Polytechnic  Research  &  Development  Co.  Inc.  1 
202  Tillary  Street 
Brooklyn  1,  New  York 
Attns  Technical  Library- 

Radiation  Engineering  Laboratory  1 

Main  Street 

Maynard,  Massachusetts 

Attns  Dr.  John  Ruze 

Radiation,  Inc.  1 

Po  0.  Drawer  37 
Melbourne,  Florida 

Attn:  Technical  Library,  Mr.  M.  L.  Cox 

Radiation  Systems,  Inc.  1 

I4I4.O  Swann  Avenue 
Alexandria,  Virginia 
Attn:  Library 

RCA  Laboratories  1 

David  Sarnoff  Research  Center 

Princeton,  New  Jersey 

Attn:  Miss  Fern  Cloak,  Librarian 

Radio  Corporation  of  America  1 

Defense  Electronic  Products 
Building  10,  Floor  7 
Camden  2,  New  .Jersey 

Attns  Mr.  Harold  J.  Schrader,  Staff  Engineer 
Organization  of  Chief  Technical 
Administrator 

Radio  Corporation  of  America  1 

Surface  Communications  Systems  Laboratory 
75  Varick  Street 
New  York  13,  New  York 
Attn:  Mr.  S.  Krevsky 

Radio  Corporation  of  America  1 

West  Coast  Missile  and  Surface  Radar  Division 

Engineering  Library,  Building  306/2 

Attns  L„  R.  Hund,  Librarian 

8500  Balboa  Boulevard 

Van  Nuys,  California 
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Code 

I  930 

I  292 

I  9x1 

I  373 

I  29h 

I  k!2 

I  295 

I  11*2 


Organization 


Radio  Corporation  of  America 
Defense  Electronic  Products 
Advanced  Military  Systems 
Princeton,  New  Jersey 
Attn:  Mr,  David  Shore 

Director,  USAF  Project  RAND 
Via:  AP  Liaison  Office 
The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 

The  Rand  Corporation 
1700  Main  Street 
Santa  Monica,  California 
Attn:  Technical  Library 

Rantec  Corporation 

23999  Ventura  Boulevard 

Calabasas,  California 

Attn:  Grace  Keener,  Office  Manager 

Raytheon  Company 

Wayland  Laboratory 

Wayland,  Massachusetts 

Attn:  Miss  Alice  G.  Anderson,  Librarian 

Raytheon  Company 
Missile  Systems  Division 
Hartwell  Road 
Bedford,  Massachusetts 
Attn:  Donald  Ho  Archer 

Republic  Aviation  Corporation 
Farmingdale,  Long  Island,  New  York 
Attn:  Engineering  Library 

Sanders  Associates,  Inc, 

95  Canal  Street 
Nashua,  New  Hampshire 
Attn:  Mr,  Norman  R.  Wild 


No,  of  Copies 

1 

1 

1 

1 

1 

1 

1 

1 
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Code 

Organization 

No. 

of  Copies 

I  96 

Sandia  Corporation,  Sandia  Base 

P«  Oo  Box  5800 

Albuquerque,  New  Mexico 

Attn:  Classified  Document  Division 

1 

I  546 

Sandia  Corporation 

Attn:  Organization  1423 

Sandia  Base 

Albuquerque,  New  Mexico 

1 

I  682 

Scanwell  Laboratories,  Inc. 

6601  Scanwell  Lane 

Springfield,  Virginia 

1 

I  312 

Space  Technology  Laboratories,  Inc. 

P0  0„  Box  95001  _ 

Los  Angeles  45.,  California 

Attn:  Technical  Information  Center 
Document  Procurement 

1 

I  297 

Sperry  Gyroscope  Company 

Great  Neck,  Long  Island,  New  York 

Attn:  Florence  ¥.  Turnbull 

Engineering  Librarian 

1 

I  367 

Stanford  Research  Institute 

Documents  Center 

Menlo  Park,  California 

Attn:  Acquisitions 

1 

I  10U 

Syl vania  Electric  Products,  Inc. 

100  First  Avenue 

Waltham  54 ,  Massachustts 

Attn:  Charles  A„  Thornhill,  Report  Librarian 
Waltham  Laboratories  Library 

1 

I  818 

Sylvania  Reconnaissance  Systems  Lab. 

Box  188 

Mountain  View,  California 

Attn:  Marvin  D0  Waldman 

1 

AF-4561/H 
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Code 
I  2 1*0 

I  338 

I  708 

I  139 

U  61 

U  15>7 

U  99 

U  3 


-  Page  17 


Organization 


TRG,  Inc0 

2  Aerial  Way 

Syosset,  New  York 

Attn:  M»  L.  Henderson,  Librarian 

A,  S.  Thomas ,  Inc, 

355  Providence  Highway 
Westwood,  Massachusetts 
Attn:  A,  S.  Thomas,  President 

Texas  Instruments,  Inc , 

6000  Lemmon  Avenue 
Dallas  9,  Texas 
Attn:  John  B.  Travis 

Systems  Planning  Branch 

Westinghouse  Electric  Corp. 

Electronics  Division 
Friendship  Int'l  Airport  Box  1897 
Baltimore  3,  Maryland 
Attn:  Engineering  Library 

Brown  University 

Department  of  Electrical  Engineering 
Providence,  Rhode  Island 
Attn:  Dr,  C„  M,  Angulo 

California  Institute  of  Technology 
Jet  Propulsion  Laboratory 
1*800  Oak  Grove  Drive 
Pasadena,  California 
Attn:  Mr,  I„  E.  Newlan 

California  Institute  of  Technology 
1201  E,  California  Street 
Pasadena, ' California 
Attn:  Dr,  C.  Papas 

Space  Sciences  Laboratory 
Leuschner  Observatory 
University  of  California 
Berkeley  1*,  California 
Attn:  Dr,  Samuel  Silver,  Professor  of 
Engineering  Science  and  Director, 
Space  Sciences  Laboratory 


No,  of  Copies 

1 

1 

1 

1 

1 

1 

1 

1 
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Code 

U  100 

II  289 

U  239 

U  183 

U  238 

U  10 

U  86 


Organization 


University  of  California 
Electronics  Research  Lab. 

332  Cory  Hall 
Berkeley  k3  California 
Attn:  J.  R.  Whinnery 

University  of  Southern  California 
University  Park 
Los  Angeles,  California 
Attn:  Dr.  Raymond  L.  Chuan 

Director,  Engineering  Center 

Case  Institute  of  Technology 
Electrical  Engineering  Department 
10900  Euclid  Avenue 
Cleveland,  Ohio 

Attn:  Professor  Robert  Plonsey 

Columbia  University 
Department  of  Electrical  Engineering 
Morningside  Heights,  New  York,  New  York 
Attn:  Dr.  Schlesinger 

University  of  Southern  California 
University  Park 
Los  Angeles  7,  California 
Attn:  Z.  A.  Kaprielian 

Associate  Professor  of  Electrical 
Engineering 

Cornell  University 

School  of  Electrical  Engineering 

Ithaca,  New  York 

Attn:  Professor  G.  C.  Dalman 

University  of  Florida 

Department  of  Electrical  Engineering 

Gainesville,  Florida 

Attn:  Professor  M.  H.  Latour,  Library 


No.  of  Copies 

1 

1 

1 

1 

1 

1 

1 
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Code 
U  59 

U  102 

II  103 

U  10k 

U  169 

U  2li0 

U  22 


Organization 


Library 

Georgia  Technology  Research  Institute 
Engineering  Experiment  Station 
722  Cherry  Street*  N.W. 

Atlanta,  Georgia 

Attn:  Mrso  J«  H«  Crosland*  Librarian 

Harvard  University 
Technical  Reports  Collection 
Gordon  McKay  Library 
303  Pierce  Hall 
Oxford  Street 

Cambridge  38*  Massachusetts 
Attn:  Librarian 

University  of  Illinois 
Documents  Division  Library 
Urbana*  Illinois 

University  of  Illinois 
College  of  Engineering 
Urbana*  Illinois 

Attn:  Dr„  P„  E„  Mayes*  Department  of 
Electrical  Engineering 

Illinois  Institute  of  Technology 
3301  S„  Dearborn  Street 
Chicago  16*  Illinois 
Attn:  Dr„  George  I.  Cohn 

Illinois  Institute  of  Technology 
Technology  Center 

Department  of  Electrical  Engineering 
Chicago  l6*  Illinois 
Attn:  Paul  C„  Yuen 

Electronics  Research  Laboratory 

The  John  Hopkins  University 
Homewood  Campus 
Baltimore  18*  Maryland 
Attn:  Dr,  Donald  E.  Kerr 

Department  of  Physics 


No«  of  Copies 

1 

1 

1 

1 

1 

1 
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Code 

U  10$ 


U  228 


U  68 


U  32 


U  26 
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The  John  Hopkins  University  1 

Applied  Physics  Laboratory 

8621  Georgia  Avenue 

Silver  Spring,  Maryland 

Attn:  Mr.  George  L.  Seielstad 

University  of  Kansas  1 

Electrical  Engineering  Department 
Lawrence,  Kansas 
Attn:  Dr.  H.  Unz 

Lowell  Technological  Institute  1 

Research  Foundation 

P.  0.  Box  709 

Lowell,  Massachusetts 

Attn:  Dr.  Charles  R.  Mingins 

Massachusetts  Institute  of  Technology  1 

Research  Laboratory  of  Electronics 

Building  26,  Room  327 

Cambridge  39,  Massachusetts 

Attn:  John  H.  Hewitt 

Massachusetts  Institute  of  Technology  1 

Lincoln  Laboratory 
P.  0.  Box  73 

Lexington  73,  Massachusetts 
Attn:  Mary  A.  Granese,  Librarian 

McGill  University  1 

Montreal,  Canada 

Attn:  Professor  G.  A.  VJoonton 

Director,  The  Eaton  Electronics 
Research  Laboratory 

University  of  Michigan  1 

Electronic  Defense  Group 

Engineering  Research  Institute 

Ann  Arbor,  Michigan 

Attn:  J.  A.  B0yd,  Supervisor 

University  of  Michigan  1 

Engineering  Research  Institute 

Willow  Run  Laboratories,  Willow  Run  Airport 

Ypsilanti,  Michigan 

Attn:  Librarian 
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U  108 


U  19k 


U  39 


U  96 


U  78 


U  109 


U  1*5 

U  185 


University  of  Minnesota  1 

Minneapolis  l4>  Minnesota 

Attn:  Mr,  Robert  H.  Stumm,  Library 

Physical  Science  Laboratory  1 

New  Mexico  College  of  Agriculture 

and  Mechanic  Arts 

State  College,  New  Mexico 

Attn:  Mr,  H.  W,  Haas 

New  York  University  1 

Institute  of  Mathematical  Sciences 

Room  802,  25  Waverly  Place 

New  York  3,  New  York 

Attn:  Professor  Morris  Kline 

Northwestern  University  1 

Microwave  Laboratories 
Evanston,  Illinois 
Attn:  R.  E,  Beam 

Ohio  State  University  Research  Foundation  1 

13lU  Kinnear  Road 
Columbus  8,  Ohio 
Attn:  Dr,  T.  E.  Tice 

Department  of  Electrical  Engineering 

The  University  of  Oklahoma  1 

Research  Institute 
Norman,  Oklahoma 

Attn:  Professor  C,  L.  Farrar,  Chairman 
Electrical  Engineering 

The  Pennsylvania  State  University  1 

Department  of  Electrical  Engineering 
University  Park,  Pennsylvania 

University  of  Pennsylvania  1 

Institute  of  Cooperative  Research 

3U00  Walnut  Street 

Philadelphia,  Pennsylvania 

Attn:  Department  of  Electrical  Engineering 
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U  48 


U  97 


U  184 


U  176 


U  110 


u  309 


U  186 


U  111 


Polytechnic  Institute  of  Brooklyn  1 

Microwave  Research  Institute 

55  Johnson  Street 

Brooklyn*  New  York 

Attn:  Dr.  Arthur  A.  Oliner 

Polytechnic  Institute  of  Brooklyn  1 

Microwave  Research  Institute 

55  Johnson  Street 

Brooklyn*  New  York 

Attn:  Mr.  A.  E,  Laemmel 

Purdue  University  1 

Department  of  Elec trie  si  Engineering 
Lafayette*  Indiana 
Attn:  Dr.  Schultz 

Stanford  University  1 

¥„  ¥„  Hasen  Laboratory  of  Physics 
Stanford*  Califo  rnia 
Attn:  Microwave  Library 

Syracuse  University  Research  Institute  1 

Collendale  Carpus 
Syracuse  10*  New  York 
Attn:  Dr.  C„  S.  Grove*  Jr. 

Director  of  Engineering  Research 

Technical  University  1 

Oestervoldgade  10  G 

Copenhagen*  Denmark 

Attn:  Professor  Hans  Lottrup  Knudsen 

University  of  Tennessee  1 

Ferris  Hall 
¥0  Cumberland  Avenue 
Knoxville  16*  Tennessee 

The  University  of  Texas  1 

Electrical  Engineering  Research  Lab. 

P.  0.  Box  8026 
University  Station 
Atustin  12*  Texas 

Attn:  Mr,  John  R.  Gerhardt,  Assistant  Diector 
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U  51 

U  132 

U  133 

U  187 


Organization 


The  University  of  Texas 
Defense  Research  Laboratory 
Austin,  Texas 

Attn:  Claude  W.  Horton,  Physics  Library 

University  of  Toronto 
Department  of  Electrical  Engr. 

Toronto,  Canada 

Attn:  Professor  G.  Sinclair 

University  of  Washington 
Department  of  Electrical  Engineering 
Seattle  5*  Washington 
Attn:  G.  Held,  Associate  Professor 

University  of  Wisconsin 
Department  of  Electrical  Engineering 
Madison,  Wisconsin 
Attn:  Dr.  Scheibe 
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